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Abstract 

Muon-spin relaxation (/x"'"SR) is a sensitive probe of magnetism, but its utility can be severely 
limited by the lack of knowledge of the muon implantation site and the extent to which the 
muon perturbs its host. We demonstrate systematically that these problems can be addressed 
accurately using electronic-structure calculations. We show that diamagnetic muons introduce 
significant short-ranged distortions in ionic insulators that would lead to systematic errors on 
magnetic moments determined by //"''SR, and quantify these. The bound state formed 

by muons in many fluorides can be understood as an exotic molecule-in-a-crystal defect with a 
zero-point energy larger than that of any naturally-occurring triatomic molecule. 
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Muon-spin relaxation (/x+SR) involves implanting spin-polarized positive muons in a sam- 
ple in order to probe the local magnetic structure [T]. /i"'"SR is an extremely sensitive probe 
of magnetism [2] but has two significant limitations. The first concerns the lack of knowl- 
edge of the site of implantation of the muon, which hinders the measurement of magnetic 
moments using /i+SR. Second, the unknown extent of the perturbation due to the muon on 
the local crystal and electronic structure of the host has been the cause for increased concern 
since /x"'"SR is frequently employed in the study of systems that lie on the verge of order- 
ing [2H3] or where doping is a critical parameter [BHH] . Previous first- principles studies have 
focussed on the paramagnetic states formed by muons and protons in semiconductors P^TT] . 
Diamagnetic muon states (where the contact hyperfine coupling is negligible) have received 
considerably less attention, in spite of their greater utility in the study of magnetic materi- 
als. Here we present a detailed microscopic study based on density-functional theory (DFT) 
of the dia- and paramagnetic muon states in a series of fiuorides, where detailed information 
about the geometry of the diamagnetic muon site is experimentally accessible, enabling an 
accurate comparison with first-principles predictions. 

In host compounds containing fiuorine, diamagnetic muons can couple strongly to the 
fiuoride ions often forming linear F-fi-F bound states \l2i, although bent F-/i-F and F-/x 
geometries have been shown to exist as well [13j. The magnetic dipolar coupling between 
muon and fiuorine nuclear spin gives rise to a signal that is sensitive to the geometry of the 
muon-fiuorine state allowing an accurate experimental determination of the muon's local 
site geometry [121 US]- In the series of non-magnetic ionic insulators LiF and NaF (rock- 
salt structure, a = 4.03 A and 4.78 A), CaF2 and BaF2 (fiuorite structure, a = 5.46 A and 
6.2 A), and for the antiferromagnetic (AFM) insulator C0F2 (rutile-type structure, a = 4.7 A 
and c = 3.18 A), we demonstrate the high accuracy with which quantitative information 
about the muon both in the dia- and the paramagnetic state can be obtained with DFT. 
We show that diamagnetic muon states lead to significant short-ranged perturbations of the 
host, not limited to the fiuorides bound in the F-/i-F states (indeed the distortion of the 
neighbouring cations can exceed that of the fiuorides). This introduces systematic errors 
on magnetic moments in ionic insulators determined by yU^SR, which we quantify from first 
principles. We study the quantum behaviour of the muon and the much heavier proton in 
both charge states. 

The ab initio calculations were performed with the Quantum ESPRESSO package [H] 
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within the generahzed-gradient approximation [15]. Unless indicated otherwise, calculations 
were performed in a supercell containing 2x2x2 conventional unit cells. The charge state 
of the muon was determined by the charge of the supercell (+1 for diamagnetic and neutral 
for paramagnetic states). A muon was placed in several randomly chosen low-symmetry 
sites and all ions were allowed to relax until the forces on all ions and the energy change 
had fallen below a convergence threshold [16] . 



(a) in LiF/NaF (b) Muonium in LiF/NaF 





FIG. 1: (color online). Calculated equilibrium geometries of dia- and paramagnetic muon states in 
LiF/NaF (Li/Na blue, F green), CaF2/BaF2 (Ca/Ba red), and C0F2 (Co magenta). Translucent 
spheres represent the equilibrium ionic positions before the muon (brown) is introduced into the 
crystal. Black lines are a guide to the eye. The c axis is vertical. 

We first discuss the equilibrium geometries of the muon states obtained from our cal- 
culations, shown in Fig. [T] Our results predict the formation of diamagnetic linear F-/i-F 
states in all of the compounds of this series. This is in agreement with previous experimen- 
tal data for non-magnetic LiF, NaF, CaF2, and BaF2, and the calculated F-fi-F states lie 
along the (110) set of directions as previously suggested [12]. The calculated and experi- 
mentally measured bond lengths are tabulated in table |T| Our calculated bond lengths are 
in all cases within 3% of the experimental values demonstrating the high level of accuracy 
of these results. We have found no evidence for any other stable diamagnetic state in this 
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series. For C0F2 a detailed experimental study [T7] has determined the muon site to be the 
octahedral ^00 site, in agreement with our calculations, although no experimental reports 
of an F-/i-F state exist. In that experiment [T71 dB], the paramagnetic region, where a 
potential F-^-F signal would be present, was studied only under an applied magnetic field, 
impeding the observation of a potential F-fi-F signal. Following our first-principles results, 
we have therefore searched for experimental signatures of an F-/x-F state in C0F2 in zero 
applied field. We have found unambiguous experimental evidence for a symmetric, linear 
F-/i-F state with a fiuoride-fluoride separation of 2.43(2) A (TU] in good agreement with 
our calculated value of 2.36 A. These results demonstrate that DFT is a powerful tool for 
determining diamagnetic muon sites. 

Using a supercell approach, we can also quantify the extent of the perturbation of the 
implanted muon on its host. The calculated structures (Fig. [T]) allow us to study the radial 
displacement of the ions away from the muon as a function of their unperturbed distance 
from the muon site (Fig. [2]). The calculated displacements demonstrate that the muon's 
perturbation is large but short ranged. While it is known that the perturbation of the 
fiuoride ions must be significant based on the experimentally measured F-/i bond lengths of 
the F-/i-F states found in many fiuorides [121 [13] , we can now quantify the perturbation of 
the cations as well. Since localized magnetic moments would be located on the cation, the 
cation displacements are particularly pertinent to understanding the effect of the muon's 
perturbation on experimentally measured //"""SR spectra discussed below. Our results show 
that in LiF, CaF2, and BaF2 the perturbation of the nearest neighbour (n.n.) cations 
even exceeds those of the fluoride ions bound in the F-/i-F state. At short distances the 
direct Coulomb interaction between the muon and the surrounding ions dominates over the 
elastic interaction transmitted through the lattice. We therefore expect similar distortions 
to be present in any ionic insulator, regardless of whether it contains fluoride ions or not. 
Indeed we believe that the formation of the F-/i-F state somewhat mitigates the n.n. cation 
distortions due to the attraction of negative charge density towards the muon. At short 
distances all displacements are radially in- or outwards from the muon due to the symmetry 
of the site. Beyond the n.n. shell, elastic interactions cause some non-radial displacement. 

We have also investigated the effect that the muon has on the magnetic moments of 
the surrounding ions. In this series only C0F2 is magnetic. The spin-only moment was 
estimated from a Lowdin population analysis to be 2.68/iB per Co ion. This compares with 
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FIG. 2: (color online). The radial displacement of the ions from the defect against the distance of 
the ion from the muon site in an unperturbed crystal. Black circle (cation), blue square (F). For 
all compounds the displacements are well converged already on a 2 x 2 x 2 supercell but for LiF 
and C0F2 the F-//-F displacements are shown for larger cells. 

a total moment of 2.60(4)/iB measured with powder neutron diffraction [2T] and a spin-only 
moment of 2.21(2)/iB determined from high energy photon diffraction [22]. We find the 
largest perturbation of the Co spin-only moment due to the presence of the diamagnetic 
muon to be about ±0.5% and therefore negligible. We believe that the perturbation of the 
total moment will be similar and therefore have a negligible effect on experimental /^"'"SR 
spectra 

In a yU"'"SR experiment, a muon at position couples to the magnetic moments of the host 
via the dipolar interaction B2[p{r^) = D'^^ {r f^)m^ , where nij is the magnetic moment of 
ion i (assumed to be completely localized) at position Vi and = k-2 ' — 5"") is the 

i i 
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TABLE I: Calculated (DFT) and experimental (exp) properties of the diamagnetic F-^-F states 
in solid and vacuum, and of the (FHF)~ molecular ion in vacuum, r (A) is the muon- fluoride bond 
length, u is the frequency (cm^-*^) of the symmetric stretch (SS), asymmetric stretch (AS), and 
bending (B) mode, and ZPE is the zero-point energy (eV). a: Our calculation, b: experimental 
data from Ref. [TO], c: Ref. [50] reports 1377 cm^^, d: from Ref. [12] 

dipolar tensor with Rj = (i?^', i?f , i?^) = r^j — rj. There is negligible contact hyperfine cou- 
pling for the diamagnetic F-/i-F muons, so in C0F2 the F-/i-F muons only probe i?jip(r^), 
which we have calculated for both a perfect crystal (i.e. in the absence of a muon) and in an 
aperiodic perturbed crystal (simulating the presence of the muon) [16j. The perturbation 
of the magnetic moments has been neglected. Our calculations predict a reduction of the 
dipolar field at the muon site by 21.4% (2x2x2 supercell), 23.6% (3x2x3 supercell), 
23.9% (4x2x2 supercell). Experimentally the dipolar coupling is measured to be 16% 
lower [17J than expected from a Co moment of 2.64/iB. This is in reasonable agreement with 
our prediction and demonstrates that relaxed geometries obtained from DFT are suitable 
for calculating corrections to expected dipolar fields and hence magnetic moments measured 
using /i+SR. Note that the perturbation of the n.n. cations in C0F2 is fairly moderate in 
comparison with the other compounds in this series (Fig. [2]). Nonetheless the n.n. cation 
displacements in C0F2 have a significant effect on the calculated dipolar coupling due to the 
short-ranged nature of the dipolar interaction [16j. This illustrates that in ionic insulators 
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the muon's perturbation on its host cannot be neglected if magnetic moments are to be 
measured accurately by yU+SR. 

All of the results above are valid for a positive muon, a proton, or a deuteron defect 
and are 'classical' in the sense that they do not take account of zero-point effects due to 
the small muon mass (m^ ~ mp/9). We have used density- functional perturbation theory 
(DFPT) [23j to calculate the vibrational properties of the F-fi-F system in the solid and 
in vacuum (table |l]). In vacuum the linear (F-/i-F)~ anion has four vibrational modes: 
symmetric stretch, bending (two-fold degenerate), and asymmetric stretch. In the solid the 
two-fold degeneracy of the bending mode is broken due to the symmetry of the site; in LiF, 
NaF, CaF2, and BaF2 there are two fundamentally inequivalent directions of bending, one 
is towards a neighbouring cation and is shifted up in frequency while the other direction 
is into a 'gap' in the crystal structure and is shifted down in frequency. The asymmetric 
stretch is very similar to its vacuum value except in CaF2, where the small bond length 
leads to a larger value. In CaF2, BaF2, and C0F2 all F-/i-F modes are highly-localized 
and decouple from the lattice modes. The decoupling of the vibrational modes of the linear 
F-fi-F system illustrates that it may be viewed as a molecule-in-a crystal defect similar to 
the Vk center found in the alkali halides [21]. In LiF and NaF the symmetric stretch mode 
mixes with the lattice modes so this analogy is slightly less apposite. From the frequencies 
of the decoupled vibrational modes we have estimated the zero-point energy (ZPE) of the 
system (in the harmonic approximation). Due to the strong muon-fluoride bond and the 
small mass of the muon, the F-fi-F center possesses an exceptionally high ZPE of 0.80 eV 
in vacuum, larger than that of any natural triatomic molecule (the ZPEs of H2O and 
are 0.56 and 0.54 eV respectively [25]). 

We now discuss the properties of the neutral muonium (Mu) state in this series. The 
calculated equilibrium geometries are shown in Fig. [1} Since the muon charge is screened by 
an electron, the Mu site is fundamentally different to the diamagnetic site. In LiF and NaF, 
Mu occupies the octahedral interstitial site. An interstitial site was previously suggested 
based on the observed hyperfine coupling [26] . In CaF2 and BaF2, Mu occupies the octahe- 
dral cation-cation centered site, the fluoride-fluoride centered octahedral site is unstable. In 
C0F2 we find a single Mu site in a nearly face-centered position at approximately (0.56, 0.84, 
0.50), distorted by the neighbouring fluoride. Fig. |2] shows the radial displacements of the 
ions from the Mu defect. Due to the screening effect of the Mu electron, the displacements 
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are generally much lower than for the diamagnetic /i"*". In LiF, NaF, CaF2, and BaF2 the 
displacements of the n.n. ions are again along the radial direction due to the symmetry of 
the Mu site. In these compounds the ions are only displaced away from the Mu. In C0F2 the 
symmetry of the Mu site is lower leading to small displacements in the tangential direction, 
even for the n.n. shell, through elastic interactions with the lattice. These are also the likely 
cause for the effective attraction of some ions despite the neutral charge state of Mu. Note 
that all of this also applies to neutral interstitial hydrogen Hp 
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TABLE II: Calculated and experimental contact hyperfine couplings (MHz) and zero-point energies 
E (eV) of Mu and H defects. A: 'classical' hyperfine coupling; {A): quantum corrected value using 
the harmonic approximation (HA), by solving the full Schrodinger equation using finite differences 
(FD). The ratio of -Era for the Mu and H modes is close to \Jmp/m^ ~ 3 indicating highly- 
localized modes, a: From Ref. |26j, b: at 295 K 127], c: at 38 K d: at 10 K [22], e: measured 
to be 0.12 eV at 100 K [30], f: at 77 K [3l], g: at 295 K [32], h: while no experimental value ^exp 
has been reported for C0F2, a total (dipolar+contact) coupling of 1280 MHz has been measured in 
the related MnF2 [33] • The dipolar coupling (without quantum correction) in C0F2 is « 71 MHz 
(0.52 T along c) and adds to the contact term. 
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The paramagnetic state is experimentally characterized by the (dipolar and contact) hy- 
perfine coupling between the muon (proton) spin and the surrounding spin density. For 
all paramagnetic states above, except the one in C0F2, the dipolar coupling cancels by 
symmetry. The dipolar coupling with the Co moments was estimated as above. Unlike in 
the diamagnetic case, the n.n. Co spin-only moment is significantly perturbed (—25%) by 
the presence of the Mu. The estimate of the dipolar coupling assumes that the relative 
perturbation of the total moment is also —25% and takes account of crystallographic dis- 
tortions [in]. The contact hyperfine coupling A is related to the unpaired spin density p(r„) 
at the muon/proton position r„ via A = ^7e7nP(rn); where 7e is the electron gyromag- 
netic ratio and 7„ is the muon/proton gyromagnetic ratio. The spin density was obtained 
using the projector- augmented wave reconstruction method j3l] and the resulting contact 
hyperfine couplings are shown in table [11} However, due to the ZPE of the defect, the defect 
wavefunction has a finite spread leading to a quantum correction to the hyperfine coupling. 
This has been previously studied in Si, Ge, and diamond (TUl [H]. Although a complete 
treatment would involve a parameterization of the full three-dimensional contact hyperfine 
coupling and potential energy surface to calculate the three-dimensional wavefunction |T0] . 
here we obtain an estimate of this correction as follows. The vibrational modes of the defect 
were calculated from DFPT and the potential energy and hyperfine coupling were calculated 
along the direction of the eigenmodes. Since the eigenmodes are mutually perpendicular, 
the motion along the three modes decouples and the wavefunction factorizes. We have then 
calculated the Mu and H? wavefunctions in two ways. The first is an anisotropic harmonic 
approximation where the wavefunction along each mode is the ground state wavefunction of 
the harmonic oscillator with the frequency u given by the calculated vibrational frequency. 
Inspection of the potential energy surface has revealed significant anharmonic terms along 
some directions in some of the compounds [T6]. We have therefore also solved the full 
Schrodinger equation along the calculated mode directions using a finite differences method. 
From the calculated contact hyperfine couplings A{r) and the Mu/H? wavefunction ip{r), we 

r dr\ip(r)\'^ A(r) 

have obtained the estimated quantum correction (A) from (A) = J r^dr\ip{r)\^ — ' "^^^ calcu- 
lated couplings A{r) and wavefunction ip{r) are weighted by to obtain an approximate 
three-dimensional average. This approximation is accurate so long as ipi^) ^ind A{r) are 
approximately spherically symmetric between neighbouring points on the integration grid. 
Since the average is taken over three mutually perpendicular directions, we expect this to be 
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a reasonable approximation. The quantum corrections are also tabulated in table [TTj Note 
that this correction is smaller for the heavier H?. Our DFT estimates are within 10% of the 
experimental value for LiF and NaF, the same level of accuracy as previous calculation in 
Si, Ge, and diamond [101 E]) ^ind within 2% of the experimental value for CaF2. 

In summary we have demonstrated systematically how DFT can be used to comprehen- 
sively address the two most fundamental limitations of the /i"'"SR technique: the problem of 
the unknown muon site and the perturbation exerted by the muon on its host. We note that 
the detailed understanding of the nature of the muon's state in solids is relevant beyond 
the field of yU+SR since the muon acts as a light analogue of hydrogen, which is a ubiqui- 
tous impurity in all technologically important semiconductors, where it strongly affects the 
electronic and structural properties of the material |35j . 
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Supplemental Information 
A. Computational details 

The ab initio calculations were performed with the Quantum ESPRESSO package |Slj 
within the generalized-gradient approximation (GGA) |S2j using norm-conserving and ultra- 
soft [S3j pseudopotentials. The functional did not include the spin-orbit interaction. The 
muon was modelled by a norm-conserving hydrogen pseudopotential. Total energies were 
converged to at least 2 x 10^^^ Ry/atom (Ry is the Rydberg constant). Structural relax- 
ations, total energies, band structures, and vibrational modes were calculated with wave- 
function and charge density cutoffs of 80 and 320 Ry. Contact hyperfine couplings were 
calculated using the projector- augmented wave (PAW) method [S4] as implemented in the 
GIPAW package |Slj . The PAW calculation required norm-conserving datasets with wave- 
function and charge density cutoff of 120 Ry and 480 Ry, respectively. The PAW datasets 
and cutoffs were also used for the calculation of the Lowdin charges used to estimate mag- 
netic moments. We used a Gaussian broadening of the occupations of 0.01 Ry. Vacuum 
vibrational modes and hyperfine couplings were calculated in a 20 x 20 x 20 A'^ cell. The 
convergence thresholds for the structural relaxation were: 10^'^ Ry/ao for the forces (oq is 
the Bohr radius) and 10"'^ Ry for the energy change between subsequent steps. The calcu- 
lated lattice parameters for the bulk compounds are within 2% of the experimental values. 
The experimental values were used for subsequent calculations. The unit cell was not varied 
during the structural relaxations which included the muon. By symmetry the ions cannot 
relax from their experimental positions in the bulk for LiF, NaF, CaF2 and BaF2. In C0F2 
we found a small relaxation of the fluorides by 0.037 A with respect to their experimental 
positions. Figs. 1 and 2 in the main text use the calculated fluoride positions for the bulk 
as reference. 

The considerable reduction of the spin-only moment in C0F2 from the ideal S = 3/2 value 
of 3/iB has been related |S5j to the relative size of the E and D anisotropy constants, where 
the anisotropy term is of the form DSl + E{Sl — Sy). Crystalline anisotropics are not repro- 
duced accurately without explicit inclusion of spin-orbit coupling in the GGA functional. 
This is the likely cause for the overestimate of the spin-only moment in our calculations 
[calculated: 2.68/xb per Co ion, high-energy photon diffraction [HE]: 2.21(2)/xb]- A previous 
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DFT study using the GGA |S7] has found a similar spin-only moment of 2.62/iB. Our cal- 
culations also predict a spin-only moment of 0.026/xb on each fluoride ion that experiences 
a perturbation of up to ±0.017;Ub due to the muon. Since the fluorine moments are very 
small, they have been neglected in the calculation of the dipolar coupling. 

In the calculation of the dipolar coupling in the aperiodic perturbed crystal, both for the 
diamagnetic and the paramagnetic state, the perturbation of the crystal structure due to 
the muon was included within one supercell around the muon (with the muon at the center), 
surrounded by unperturbed ions within a radius of 23 a lattice parameters. 

In general, the magnetic coupling of the muon is a sum of dipolar coupling, contact 
hyperfine interaction, demagnetization and Lorentz fields. The calculated contact hyperfine 
coupling for the diamagnetic F-/i-F muon in C0F2 is < 1 MHz and is therefore negligible in 
comparison with the dipolar coupling. In an antiferromagnet the demagnetizing and Lorentz 
fields are zero. Hence for the diamagnetic F-/i-F muon only the dipolar coupling needs to be 
considered. The short-ranged nature of the dipolar coupling can be illustrated by separately 
considering the dipolar coupling of the F-/i-F muon with the nearest neighbour (n.n.) and 
the remaining cations in the crystal: assuming a Co moment of 2.64/iB, the crystallographic 
distortion reduced the dipolar field from the two n.n. Co ions from —0.378 to —0.309 T, 
while the dipolar field from all other ions in the solid merely changes from 0.113 to 0.101 T 
(all along the c-axis and for the 2x2x2 supercell). If the cations in the non-magnetic 
compounds LiF, NaF, CaF2, and BaF2 were magnetic and aligned antiferromagnetically 
(ferromagnetically) along the c-axis, the dipolar corrections for the F-/i-F states shown in 
Fig. 1 in the main text would be: LiF -39%(-70%), NaF -34%(-67%), CaFs 0(-46%), and 
BaF2 0%(-47%). In the antiferromagnetic case, the muon would be located in a site of 
cancellation of the dipolar field in CaF2 and BaF2. 

For Mu in C0F2 the dipolar coupling was estimated from the dipolar coupling to the Co 
moments only. The Mu electron spin density is approximately spherically symmetric and 
therefore only yields a small contribution to the dipolar coupling, which has been neglected. 
The Co moment was assumed to be 2.64/iB and the perturbation of the n.n. moment was 
taken to be —25% (the calculated reduction of the spin-only moment). The perturbation 
of the spin-only moments of the other Co ions in the supercell was negligible. At different 
levels of approximation the dipolar coupling is 0.49 T (unperturbed crystal), 0.72 T (crys- 
tallographic distortions only) and 0.52 T (crystallographic distortions and perturbation of 
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the n.n. Co moment), all are along c and have the same sign as the contact coupling. 
B. Experimental search for F F state in cobalt (II) fluoride 
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FIG. S3: (color online). Experimental muon decay asymmetry A{t) at 70 K in C0F2. The red line 
is the fit to Eq. [Sll 



A powder sample of C0F2 (Sigma Aldrich 236128) was wrapped in 25 /im silver foil and 
mounted in a ^He cryostat on the GPS instrument at Paul Scherrer Institut in Switzerland. 
Above the critical temperature of 37.85 K, we observed oscillations in the muon decay 



asymmetry A{t) characteristic of an F-/i-F state, see Fig. S3 The data were fitted to 



A{t) = Aiexp(-Ait)L'^(t) + Abgexp(-Abgt), (SI) 

where Dz{t) describes the time evolution of the muon spin in an F-/i-F state |S8j and the 
second term accounts for a slow-relaxing background due to muons stopping in the cryostat 
tail or sample holder. The term exp(— Ait) phenomenologically takes account of residual 
magnetic dynamics in the sample with a fitted value of Ai = 0.73(8) MHz at 70 K. The best 
fit is for a symmetric, linear F-/i-F state with a fluoride-fluoride separation of 2.43(2) A in 
good agreement with our calculated value of 2.36 A. The extracted F-/i-F geometry is not 



very sensitive to the precise choice of parameterization in Eq. SI as it is dominated by the 
nature of the function Dz(t) which encodes this geometry. 
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FIG. S4: (color online). Quantum correction to contact hyperfine coupling using both the harmonic 
approximation (HA) and by solving the full Schrodinger equation using finite differences (FD) along 
two eigenmodes in C0F2. From top to bottom: schematic of the eigenmode, the red arrow indicates 
displacement in positive direction; the contact hyperfine coupling for Mu (black) and H-* (green); 
the wavefunction weighted by r^; the normalized unweighted wavefunction for Mu (HA black; FD 
blue) and H? (HA red; FD orange); the potential energy above equilibrium calculated from the 
DFT total energy and the harmonic approximation to it for Mu and H-*. While for the mode on 
the left the HA is quite reasonable, there are large deviations for the second mode. 

C. Quantum correction of contact hyperfine coupling 

In the anisotropic harmonic approximation, the ground state Mu (H?) wavefunction is 
given by 




(S2) 
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where = ri/Ri is the distance from equihbrium along eigenmode i of angular frequency uj. 
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normalized by the range of the wavefunction Ri = and m is the mass of the particle. 

Along each mode we have calculated the contact hyperfine coupling A and the total energy 
up to 3 Ri from equilibrium in positive and negative direction at equal spacings. For the 
finite differences method, we have used the total energies from the supercell calculation as 
the potential for which the Schrodinger equation was solved. The quantum average described 
in the main text is then found from 

E(W)V(rDr^(rD 

where i sums over modes and j is the step along each eigenmode. An example calculation 



is shown in Fig. S4 



D. Formation energies 

Finally we address the issue of which charge state of the muon/proton is more favourable 
at equilibrium. We define the formation energy of the defect as 

-^f = -^defect - -^Mk + -^ZPE + ^-^F, (S4) 

where -E'defectl'?) ^'^^ total energy obtained from a supercell calculation including the defect 
in charge state g, E'^^^^ is the total energy of the bulk, -Ezpe is the zero point energy of the 
defect (estimated from the decoupled vibrational modes, see main text), and is the Fermi 
energy with respect to the top of the valence band. Note that this energy scale is only relative 



and that it does not include a muon chemical potential. Fig. [85] shows the formation energy 
as a function of Fermi energy in all compounds in this series. We have used the ZPE energy 
for a muon, but the corresponding plots for hydrogen will be very similar since the differences 
in ZPE are less than 0.5 eV. In the undoped case, where the Fermi level lies halfway between 
valence and conduction band, the formation energy is lower for the neutral center in NaF 
while the diamagnetic state lies lower in CaF2, BaF2, and C0F2. LiF is approximately 
amphoteric, i.e. the muon/proton assumes a charge state that counteracts the prevailing 
conductivity. While this would explain why the neutral interstitial hydrogen center appears 
to be more favourable in NaF than in LiF |Sllt IS12j . we point out that both charge states 
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FIG. S5: (color online). Formation energies of the F-/i-F states relative to the muonium states as 
a function of Fermi energy. Calculated (-E^aic) experimental \S9\ ISIO] band gaps (E'fxp). 

have been clearly identified in /i+SR experiments, despite differences in formation energies 
of several eV. This highlights the importance of the kinematics of the stopping process for 
the formation probabilities of muon states in solids in a /i"''SR experiment in agreement with 
experimental studies using slow muons |S13j and applied electric fields |S14] . In addition, 
muons might locate in the proximity of cation vacancies not studied here. However, for 
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small vacancy concentrations, this effect is likely to be negligible. 



[SI] 
[S2] 
[S3] 
[S4] 
[S5] 
[S6] 
[S7] 
[S8] 
[S9] 
[SIO] 
[Sll] 
[S12] 
[S13] 
[S14] 



P. Giannozzi et a/., |J. Pliys.: Condens. Matter 21, 395502 (2009) 
J. P. Perdew, K. Burke, and M. Ernzerhof, |Phys. Rev. Lett. 77, 3865 (1996) 



D. Vanderbilt, jPhys. Rev. B 41, 7892 (1990) 



P. E. Blochl, iPhys. Rev. B 50, 17953 (1994) 



W. Jauch, M. Reeliuis, and A. J. Schultz, Acta Crystallogr A 60, 51 (2004) 



J. Strempfer et al, Phys. Rev. B 69, 014417 (2004) 



P. Dufek et a/., |Phys. Rev. B 49, 10170 (1994)[ 



T. Lancaster et a/.,|Phys. Rev. Lett. 99, 267601 (200"7) 



R. T. Poole, J. G. Jenkin, J. Liesegang, and R. C. G. Leckey, [P^. Rev. B 11, 5179 (1975) 
G. W. RublofF, jPhys. Rev. B 5, 662 (1972)| 



C. Hoentzsch and J. Spaeth, Sohd State Commun. 29, 577 (1979) 



T. Kamikawa, Physica Status Sohdi (b) 99, 721 (1980) 



T. Prokscha et a^., Phys. Rev. Lett. 98, 227401 (200"7)| 



V. Storchak et a/., Phys. Rev. Lett. 78, 2835 (1997) 



18 



